クロス集計
library(gmodels)
data(birthwt, package = "MASS")
CrossTable(birthwt$smoke, birthwt$low, prop.r = T, prop.c = F, prop.t = F,
prop.chisq = F, expected = T, chisq = T, fisher = T)
##
##
## Cell Contents
## |-------------------------|
## | N |
## | Expected N |
## | N / Row Total |
## |-------------------------|
##
##
## Total Observations in Table: 189
##
##
## | birthwt$low
## birthwt$smoke | 0 | 1 | Row Total |
## --------------|-----------|-----------|-----------|
## 0 | 86 | 29 | 115 |
## | 79.101 | 35.899 | |
## | 0.748 | 0.252 | 0.608 |
## --------------|-----------|-----------|-----------|
## 1 | 44 | 30 | 74 |
## | 50.899 | 23.101 | |
## | 0.595 | 0.405 | 0.392 |
## --------------|-----------|-----------|-----------|
## Column Total | 130 | 59 | 189 |
## --------------|-----------|-----------|-----------|
##
##
## Statistics for All Table Factors
##
##
## Pearson's Chi-squared test
## ------------------------------------------------------------
## Chi^2 = 4.924 d.f. = 1 p = 0.02649
##
## Pearson's Chi-squared test with Yates' continuity correction
## ------------------------------------------------------------
## Chi^2 = 4.236 d.f. = 1 p = 0.03958
##
##
## Fisher's Exact Test for Count Data
## ------------------------------------------------------------
## Sample estimate odds ratio: 2.014
##
## Alternative hypothesis: true odds ratio is not equal to 1
## p = 0.03618
## 95% confidence interval: 1.029 3.965
##
## Alternative hypothesis: true odds ratio is less than 1
## p = 0.9911
## 95% confidence interval: 0 3.581
##
## Alternative hypothesis: true odds ratio is greater than 1
## p = 0.02021
## 95% confidence interval: 1.137 Inf
##
##
##